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ABSTRACT. Semi-analytical models of galaxy formation based on hierarchi- 
cal clustering now make a wide range of predictions for observable properties 
of galaxies at low and high redshift. This article concentrates on 2 aspects: 
(1) Self-consistent modelling of dust absorption predicts a mean UV extinc- 
tion Ajjy ~ 1 mag, depending only weakly on redshift, and similar to ob- 
servational estimates. (2) The models predict that the Lyman-break galaxies 
found at z ~ 3 should be strongly clustered with a comoving correlation 
length ro = 4 — 7 Mpc/ft (depending on the cosmology), in good agreement 
with subsequent observational determinations. 



1 Introduction 

The technique of semi-analytical 
modelling is a powerful tool for mak- 
ing predictions about the observable 
properties of galaxies at low and 
high redshift, based on hierarchical 
clustering models of structure for- 
mation, and thus for testing ideas 
about galaxy formation. The models 
we describe here are a greatly im- 
proved version of those of Cole et 
al. (1994), and are fully described in 
Cole ct al. (1998). The models in- 
clude the following processes: (a) for- 
mation of dark halos through merg- 
ing; (b) shock-heating and radiative 
cooling of gas within halos; (c) col- 
lapse of cool gas to rotationally- 
supported disks; (d) star formation 
from cold gas; (e) feedback from su- 
pernovae; (f) chemical evolution; (g) 
galaxy mergers producing elliptical 
galaxies and bulges; (h) luminosity 
evolution of stellar populations; and 
(i) absorption by dust. The treat- 
ment of dust absorption is new, and 
is described in Section ^. Predictions 
for galaxy clustering at high redshift 
are described in Section tj. 

The outputs from the models are 
masses, luminosities, colours, sizes, 
morphologies, circular velocities and 
chemical compositions for the whole 
galaxy population at any redshift. 
The models depend on several pa- 
rameters relating to star formation 
and feedback, which are chosen to 
match observations of present-day 
galaxies, as described in Cole et al. 



(1998). Once normalized in this way, 
the models match most observations 
of present-day galaxies well, and can 
be used to make predictions for the 
evolution of galaxy properties with 
redshift. 



2 Dust absorption 

Absorption of starlight by dust 
has a significant effect on optical lu- 
minosities and colours of galaxies, 
and a large effect on the far-UV lumi- 
nosities which are used as the main 
tracer of star formation rates at high 
redshift. We calculate dust absorp- 
tion in a self-consistent way for each 
galaxy, using the 3D radiative trans- 
fer models of Ferrara et al. (1998). 
The latter models include stars in a 
bulge and exponential disk, and dust 
in an exponential disk, and give the 
net attenuation of galaxy light as a 
function of wavelength, inclination, 
and the galaxy parameters. These 
models are a major improvement 
over previous treatments, which as- 
sume a ID screen or slab geometry 
for the stars and dust. The most im- 
portant parameter in the dust mod- 
els is to, the central face-on optical 
depth in dust (e.g. in the V-band). 
This is calculated directly from the 
semi-analytical model, as: 

Mdust h[ ga sZgas , . 
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Figure 1. The cosmic star formation history. The curves show the model pre- 
dictions, for a CDM model with flo = 0.3, Ao = 0.7 and a Kennicutt (1983) 
IMF. The 1500A luminosity density from the model has been converted to 
a star formation density using a constant conversion factor, as for the ob- 
servations. The dashed curve shows the model without dust, and the solid 
curve the model prediction including dust. The symbols show observational 
estimates of the star formation rate per comoving volume, derived from Ha 
or UV luminosity densities assuming a Kennicutt IMF and solar metallicity. 
The triangles are from Ha luminosities, and have been corrected observa- 
tionally for dust (and so should be compared to the dashed curve), while the 
other points are from rest-frame UV (1500 — 3000^4 ) luminosities, and do not 
include any dust correction (and should be compared to the solid curve). The 
references for the observational data are: open triangle: Gallego et al. (1995) 
filled triangle: Tresse & Maddox (1998); filled square: Treyer et al. (1998) 
open circles: Lilly et al. (1996); stars: Sawicki et al. (1997); filled hexagons 
Conolly et al. (1997); open squares: Madau et al. (1998). Filled pentagons 
show preliminary results from a new survey by Steidel et al. (1998), and may 
be revised when the survey is completed. 



The dust is assumed to have a so- 
lar neighbourhood extinction law, 
with the dust /gas ratio scaling as 
the gas metallicity Z gas . Z 

gas IS 

obtained from our chemical evolu- 
tion calculation, and the disk radius 
r disk is calculated based on angu- 
lar momentum conservation. We as- 
sume a ratio of vertical scaleheights 
hz(dust)/h z (stars) = 1, but the re- 
sults are only weakly dependent on 
this value. Thus, our calculation of 
dust absorption has essentially no 
free parameters. 

We find in our models that 
the mean dust extinction increases 
strongly with galaxy luminosity, this 
being an effect mainly of the increase 
of mean surface density with galaxy 
mass. On the other hand, the de- 
pendence of net extinction on wave- 
length is much weaker than predicted 



by a simple foreground screen model 
for the dust. Figure [I] shows the ef- 
fects of dust on the inferred star for- 
mation history of the universe, for a 
CDM model with C = 0.3, A = 0.7 
and a Kennicutt (1983) IMF. Since 
most observational estimates of the 
star formation density are derived 
from far-UV luminosities, we have 
followed the same approach for the 
models, taking the predicted 1500A 
luminosity density (with or without 
dust) and converting it to an SFR 
density assuming the same IMF and 
a constant (solar) metallicity. The 
mean 1500A extinction is predicted 
to vary only slowly with redshift, in- 
creasing from 0.5 mag at z = 
to 1.2 mag at z = 6, the smaller 
dust content of higher redshift galax- 
ies being more than compensated by 
their smaller sizes, resulting in higher 



Semi-analytical galaxy formation models 



3 



optical depths. The mean extinction 
at 3000A is similar to that at 1500A. 
These values for the UV extinction 
are similar to those estimated obser- 
vationally for high redshift galaxies 
e.g. Pettini et al. (1998). 

The model shown predicts that 
the star formation rate per comoving 
volume peaks at z ~ 3. The model 
prediction appears to be broadly 
consistent with the current observa- 
tional data, but, given the large er- 
ror bars and significant scatter be- 
tween different observations, a defini- 
tive comparison is not yet possible. 



3 Clustering of high redshift 
galaxies 

Semi-analytical galaxy formation 
models predict what should be the 
galaxy populations in dark halos of 
different masses, and thus predict 
what should be the clustering bias 
of galaxies relative to dark matter. 
Measurements of the clustering of 
galaxies at different redshifts and 
different luminosities thus provide 
an important test of galaxy forma- 
tion models. The clustering ampli- 
tude has now been measured for a 
large sample of normal galaxies at 
z ~ 3, detected through their Lyman 
break features (e.g. Steidel, these 
proceedings). The predictions of our 
semi-analytical model for the prop- 
erties of these Lyman-break galax- 
ies originally found by Steidel et 
al. (1996) have been presented in 
Baugh et al. (1998). They show that 
CDM-based models with reasonable 
parameters are able to explain the 
observed number density and other 
properties of these galaxies. An im- 
portant prediction of the models, 
made in advance of the observations, 
was that the Lyman-break galaxies 
at z ~ 3 should be strongly clus- 
tered. Figure ^ shows the prediction 
of Baugh et al. (1998) for the two- 
point correlation function of Lyman- 
break galaxies brighter than appar- 
ent magnitude Rab = 25 at z fa 3, 
for two different cosmological mod- 
els. 

The Baugh et al. prediction was 
based on an analytical estimate of 
the non-linear evolution of the power 
spectrum of density fluctuations, and 
on the Mo & White (1996) analyti- 
cal formula for the clustering bias b 



of dark matter halos of mass M at 
redshift z: 



b(M, z) = 1+ 



S c (z) 



r 2 (M, Z) 



(2) 



where 8 c (z) is the critical linear over- 
density for an object to collapse at 
redshift z, and cr(M, z) is the vari- 
ance of linear density fluctuations of 
mass M. This formula for the halo 
bias is expected to be valid down to 
scales somewhat smaller than the co- 
moving radii of the halos concerned 
(shown by an arrow in Figure ^), at 
which point halo-halo exclusion ef- 
fects become significant. The semi- 
analytical models predict that the 
Lyman-break galaxies found by Stei- 
del et al. should lie in halos of masses 
M ~ 1O 12 M0. This is much larger 
than the typical halo mass M* (de- 
fined by cr(M*,z) = 5 c (z)) at z = 3, 
and so according to equation ^, these 
halos and the galaxies they contain 
should be highly biased. Baugh et 
al. calculated the mean clustering 
bias for the Rab < 25 Lyman- 
break galaxies, and found b f» 4 for 
both Q = 1 and C = 0.3, A = 
0.7 CDM models. The accuracy of 
the analytical approximations used 
for the halo clustering was subse- 
quently confirmed by Governato et 
al. (1998), who combined the semi- 
analytical models with N-body sim- 
ulations, and found essentially iden- 
tical results for the clustering of Ly- 
man break galaxies. This is shown in 
Figure |i| 

A useful measure of the clus- 
tering amplitude is the correlation 
length ro, defined to be the separa- 
tion (in comoving coordinates) where 
5( r o) = 1- The Baugh et al. pre- 
diction for the Lyman-break galax- 
ies was ro = 4 Mpc/h for Q = 
1 and ro = 7 Mpc/h for Qo = 
0.3, Ao = 0.7. These predictions 
agree remarkably well with the val- 
ues subsequently measured by Adel- 
berger et al. (1998) from a counts- 
in-cells analysis of redshift surveys 
of these galaxies, which were ro = 
(4 ± 1) Mpc/h for CI = 1 and r = 
(6±1) Mpc/h for C = 0.3, A = 0.7. 
These values are also plotted in Fig- 
ure ^ for comparison with the mod- 
els. This successful prediction of the 
clustering of high-redshift galaxies is 
strong evidence for these galaxies be- 
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Figure 2. The clustering of Lyman-break galaxies brighter than Rab = 25 at 
z as 3, based on an analytical halo clustering model. The solid line shows the 
two-point correlation function in comoving coordinates for a CDM model with 
= 1, while the dashed line shows a CDM model with Qq = 0.3, Ao = 0.7. 
The symbols with error bars show the values of ro ( defined by £ (ro) = I) found 
observationally by Adelberger et al. (1998) for different assumed cosmologies: 
filled symbol: Q = 1; open symbol: Qq = 0.3, Aq = 0.7. 




Log(r/h ^pc) 

Figure 3. The clustering of Lyman-break galaxies at z = 3, calculated using a 
combined N -body /semi- analytic method, (a) Correlation functions in comov- 
ing coordinates for the f2 = 1 CDM model. Lower solid line: dark matter cor- 
relation function; dashed line: real space correlation function of Lyman-break 
galaxies; upper solid line: redshift space correlation function of Lyman-break 
galaxies; dotted line: real space correlation function of Lyman-break galax- 
ies calculated from the analytical model, as in Figure ^. (b) Redshift space 
correlation functions for SI = I (SCDM) and O = 0.3, A = (OCDM) 
models. 
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ing in the most massive halos at this 
redshift, and is an important confir- 
mation of the validity of the semi- 
analytical models. 
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